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Abstract 

We review recent results on the behaviour of the dielectric function of cylin- 
drical nano-organic materials at very low frequencies in a magnetic field. For 
cylindrical structures - such as carbon nanotubes - the polarisability is shown 
to be a discontinuous function of a longitudinal magnetic field where plateau- 
like regions are separated by sudden jumps or peaks. A relation is pointed 
out between each discontinuity in the polarisability and the cross-over be- 
tween ground and first excited states induced by the magnetic field. This one 
to one correspondence suggests to use measurements of the dielectric func- 
tion in an applied magnetic field in order to obtain informations about the 
electronic structures of cylindrical nanostructures. In addition, it is shown, 
by studying finite graphene layers, that the measurement of the polarisability 
in a magnetic field could be a powerful way for detecting possible edge-states 
in amorphous carbon materials such as activated carbon fibres. Finally, the 
importance of the electron-electron interaction is emphasised by discussing 
examples of strongly interacting electrons on rings or cylinders, in the limit 
of infinite interaction. 
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I. INTRODUCTION 



It is known, since long time, that ring-shaped organic molecules, such as benzene or 
naphthalene, show very large diamagnetic susceptibilities and a very pronounced diamag- 
netic anisotropy: the diamagnetic susceptibility has normal values in the plane of the aro- 
matic molecules but abnormally large ones in the direction perpendicular to this plane. A 
first qualitative explanation of this phenomenon as arising from the Larmor precession of 
electrons in orbits including many nuclei can be traced back to P. Ehrenfest in 1925 M. This 
idea was further developed eleven years later by L. Pauling who proposed a quasi-classical 
description of the 7r electronic currents around the rings using an electrical network as a 
model f2|. One year latter, in 1937, F. London set up the "modern" theoretical description 
of this phenomenon by using a tight binding model which included the magnetic field in 
form of a complex phase of the usual transfer integral 0. This is the so-called London (or 
Peierls) substitution widely used nowadays. 

The comparison between the results given by these early theoretical works and the ex- 
perimental data available at that time is quite satisfactory However, the interest for 
this topic has nearly vanished most probably due to experimental limitations. For instance, 
in the theory of London, the electronic spectrum of aromatic molecules is a periodic function 
of the magnetic flux with a period equal to a flux quantum O — hc/e. But the value of 
the magnetic field needed to cover a full period is of order 10 5 T for the materials existing 
until few years ago. Such huge values cannot be reached experimentally, and are due to 
the extremely small perimeter, of order a few tens of Angstroms, of the materials studied 
previously. 

Today, material science has made considerable progress and many new carbon materials 
with nanoscopic characteristic length exist. An example is given by multi- walled carbon nan- 
otubes with a diameter of approximately 10 nm [|J. For these new materials, the magnetic 
field needed to reach a flux quantum is about a few tens of Tesla only and the periodicity 
mentioned before can be observed 0. 
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On the other hand, during the last years considerable progress has been made in precise 
measurements of the real part of the low-frequency dielectric function e'(u). In particular 
these experiments have been extended to ultra-low temperatures. Ratios of Se'/e' up to 10~ 7 
were achieved in the mK temperature regime. A rather spectacular success associated with 
that progress was the observation of a strong magnetic-field dependence of the polarisability 
of multicomponent glasses ||. These developments suggest reconsideration of the magnetic 
field effects in organic nanostructures such as large aromatic molecules or carbon nanotubes 
for instance. In particular, the magnetic field should have a very pronounced effect on the 
polarisability of these organic systems which should be observable nowadays These 
measurements could then be used as a new spectroscopic tool to obtain informations on 
the electronic properties of organic materials. The aim of this paper is to review various 
theoretical works developed in Dresden during the last two years in that direction fTf-fLOU- 

Similar ideas were already pointed out recently for quite different materials and purposes. 



Indeed, in [p^|-[T3|] small particles (rings, disks or spheres) were considered in connection with 
weak localisation. In particular, it was found that the polarisability should be greater in an 
applied magnetic field than in zero field because a field partially suppresses weak localisation. 
This theoretical prediction was observed very recently in the AC polarisability of mesoscopic 
rings 
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Although the studies already done until now on the behaviours of the dielectric suscep- 
tibility of cylindrical organic materials under magnetic field [0-0 should be improved in 
several ways by including, e.g., explicitely Coulomb interactions, we believe the phenomena 
pointed out here are robust in the sense that they will survive a more sophisticated the- 
oretical treatments. Therefore, we think that the polarisability in a magnetic field could 
be a useful quantity to study experimentally in order to obtain insight into the electronic 
structure of carbon made materials. In our opinion, this justify the topic of this review and, 
we hope that it will help to motivate new experiments. 

The series of works described here is related to the studies of persistent currents in 
metallic rings. Before starting, we would like to mention recent review articles on this topic. 



Interested readers could refer to Ref. [pl| for the case of diffusive electron motion and to 
Ref. [[16] for ballistic electron motion. 

The review is organised as follows. In section 2 we describe the Aharonov-Bohm effect in 
cylindrical nanostructures by using London's substitution. In section 3 and 4, we point out 
dramatic consequences of this effect on the dielectric function of nanoscopic organic cylinders 
such as carbon nanotubes (section 3), and a finite plane of graphite where the edge states 
existing in these materials define effective rings (section 4). Finally, before concluding, we 
point out a few properties of strongly interacting electrons on rings and cylinders in an 
applied field. 

II. AHARONOV-BOHM EFFECT IN CYLINDRICAL NANOSTRUCTURES 

In this paper, we are interested in the effects of an applied magnetic field on the orbital 
degrees of freedom only. Therefore we consider spinless fermions except in the last section. 
The Zeeman interaction was briefly discussed in |7|||. In the range of field strength we are 
interested in, the diamagnetic response is the most important one. 

Let us first consider a ring of N sites pierced by an axial magnetic field. The Hamiltonian, 
making use of the London substitution |§, is the following 

N 



n=l 



H R = tJ2(< + ia n e l ^ + h.c.) (1) 



where a£, (a n ) are spinless fermion creation (annihilation) operators and (ft is the magnetic 
flux through the ring in units of the flux quantum 0o (</>o — hc/e) . It is proportional to the 
area enclosed by the ring 

, N 2 a 2 hc TT 

where a is the lattice constant. For a ring with a few hundred sites and assuming a ~ 1AA, 
the magnetic field needed to reach <fi = 1 is not more than a few tens tesla - values reachable 
experimentally. 
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Introducing the wave function |^ >= 52%=i c n \ip n >, where \ip n > are the tc orbitals 
located at sites n, one obtains the following eigenvalue equations with periodic boundary 
conditions 

c n e = tcn+ie^ + fc n „ie-^ 

(3) 

C N+1 = Ci 

where e is the eigenvalue. By performing a pseudo-gauge transformation introduced by 
Bayers and Yang [[17]], c n = c n e l ^r^ n , Eqs. (|]) can be written in the form 



(4) 

c N+1 = Cl e l2 ^ 

The effect of the magnetic field appears in the new periodic boundary conditions. An analysis 
of them shows that the spectrum of the system is a periodic function of with period 1, 
as already mentioned in the introduction. This is an important information since it implies 
that every thermodynamic quantity is also periodic function of with the same period. 
Moreover, the system is invariant with respect to the transformation H — > —H. Therefore 
its spectrum - and every thermodynamic quantity - is symmetric with respect to <p = 0. In 
principle, it is then necessary to study the system only for e [0, 1/2]. 
The eigenvalues for the ring are 

2tt 

e p = 2tcos(— (p + (/>)) (5) 

with p — 0, ±1, ±2, .... Using this expression one notices crossing between levels occurring 
at special values of 0, namely = 0, 1/2 and 1 (for an example see Fig. 1). 

Let us now consider cylinders made by rolled square lattices with an applied axial mag- 
netic field (Fig. 2). The Hamiltonian for spinless fermions is 

N 2jr M 

Ho = t Yl (al+i, m an, m e^^ + h.c.) + t ^ (al,m+i a n,m + h.c.) (6) 

n=l m=l 

where the magnetic flux is still given by expression (§). The resulting eigenvalues are 

e M = 2t cos(^(p + 0)) + 2t cos(^yg) (7) 
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with —N/2 < p < N/2 — 1 and 1 < q < M. We have applied open boundary conditions at 
the ends of the cylinder. The resulting one-electron wave functions are 



I 9 N M 

l*» >= \/jv(mTi) „? 1 £ e '*'" sin( MTT ?m) < JO > (8) 

where |0 > is the vacuum. 

Again, using the expression (0), one notices many level crossings induced by the magnetic 
field (Fig. 3). Compared to the results obtained for rings, the important difference is that 
the crossings appear for values of depending on the geometry of the system (N,M); in 
particular, they appear for other values than 0, 1/2 and 1. Due to their intricate patterns, 
numerous changes of the ground state induced by the magnetic field do occur. As an example, 
we show in figure 4 as function of the magnetic flux the energies of the ground state and the 
first excited state of a cylinder with N = 100, M = 100 and with N e = 1000 electrons in it. 
The number of crossing points and the values of for which they occur depend crucially on 
the geometry of the system as we have already mentioned above but, also, on the electron 
filling. We will not attempt to explain the crossing point pattern here. Instead, as shown 
in the next section, we want to point out that the changes of the ground state associated 
with the crossing points, have strong influences on the static polarisability. Therefore this 
quantity could be used as a new spectroscopic tool to obtain informations on excited states 
or on the system geometry for instance. 

Carbon nanotubes form a class of cylindrical systems of growing importance nowadays 
PJT8[|. The calculations above can be straightforwardly extended to those systems, in par- 



ticular for achiral nanotubes, i.e. the armchair and zig-zag nanotubes considered here |18 



For these particular systems with a uniform magnetic field along their axis, we find by using 
again a nearest-neighbours tight binding model, the following spectrum |T8|,|7|J8[] 

e Ptq = ±(l + u p ±(u p v q ) 1/2 ) 1/2 (9) 

where for zig-zag nanotubes 
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u p = 2(1 + cos(f(p + </>))), with p=l,..,N 

«, = 2(1 + 0)8(^9)), with g=l,.,M 
and for armchair nanotubes 

m p = 2(1 + cos(-^yp)), with p=l, ..,M 

^ = 2(l + cos(f (g + 0))), with q = l,..,N 



(10) 



(11) 



The roles of u p and f g are just exchanged from one case to the other. 

Using the expressions (P- |TT|) , one notices that - as for the case of rolled square lattices 
- the ground state of rolled honeycomb lattices changes also strongly in a magnetic field. 
Therefore, as for the case of rolled-square lattices, the static electric polarisability of carbon 
nanotubes are also expected to be strongly influenced by an axial magnetic field [0|§. It 
is the purpose of the next section to analyse the structure of the dielectric susceptibility of 
cylindrical organic nanostructures in presence of an applied magnetic field. Moreover, the 
electronic spectrum of carbon nanotubes is also a periodic function of the magnetic flux, a 
behaviour which has been confirmed experimentally in |J. 

III. STATIC ELECTRIC POLARISABILITY OF CYLINDRICAL 
NANOSTRUCTURES UNDER MAGNETIC FIELD 

Very precise measurements of the polarisability can be done nowadays. Two recent 



examples are given in || for multicomponent glasses and in |14]] for mesoscopic metallic rings. 
The complicated electronic spectrum of cylindrical organic nanostructures - including carbon 
nanotubes - induced by an axial magnetic field causes important effects in the polarisability 
of these systems as we want to discuss in the following. 

For simplicity let us focus on rolled square lattices only. We give one example for carbon 
nanotubes in this section but, more informations can be found in 

With a uniform DC electric field E, the tight binding model becomes 



H = H + Ed (12) 

where d is the dipole operator. In the following, two different orientations for the electric-field 
are considered. The corresponding expression for the dipole operator are 

d = eaJ2(n~ ^-)^ n , m a n, m (13) 

n,m 

for longitudinal field, where a is the lattice spacing and L = Ma is the total length of the 
tube, and 

277 

d = eRY / cos(—n)al im a n>m (14) 

n,m 

for transversal field, where R is the radius of the cylinder (Fig. 2). The first case was studied 
in 0, the second in |§. 

We are interested in the behaviour of the induced dipole moment as function of the 
magnetic flux 

1 Trdp~P(H-Li) 

^>»4?^r < 15 > 

where, as usual, (3 = ^-^ and \i is the chemical potential which is a function of the magnetic 
field. 

We summarise in the following, the main results of the studies done in 

A. Longitudinal electric field: Semi-classical approximation 

In this case, the Hamiltonian of the system is given by Eqs. (|i~2|) and ([13]). With this 
model it becomes tedious to evaluate the polarisability (pT5|) , for long chains. Instead of 
doing exact calculations we proceed by using a semi-classical approximation very powerful 
in this specific context 0. 

We are interested in the linear regime where the potential induced by the electric field, 
eaEn, is slowly varying on the scale of the Fermi wave length. Then, the Thomas-Fermi 
approximation may be applied [jnj. We assume that it is meaningful to consider an energy 



depending on the position, m, on the chain 
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M 



~p,q 



m 



e P:q + eaE(m - — ) (16) 



where e Piq is given by the expression ([/]). Next, within this semi-classical approximation, the 
polarisability is evaluated as follows. 

P(T rh) — V ( m ~ f ) ( -x 7 \ 
{ ,<P) EM ^ e^'wH-c) + 1 1 > 

p,q,m ^ 1 

At this stage, it should be noticed that the interaction which leads to an energy splitting 
at crossing points when the electric field is turned on, is neglected here. Therefore, a possible 
first order Stark effect near those points is not taken into account with this approximation. 

To illustrate the relevance of this approximation, we consider first a linear chain in the 
absence of a magnetic field. In this case, the quasi-classical energy at position m is given by 

7T M 

e q (m) = 2tcos(^— -q) + eaE(m - —) (18) 
and the corresponding polarisability by 

P m- ea v ( m ~f) nq] 

lin{ ' EM ^ e ^( m )^) + 1 K ' 

Figure 5 shows an example for a linear chain with 201 sites. The exact polarisability 
is compared with the semi-classical one (Eq. [19]) as function of the electron filling for two 
different values of the temperature. One can see that the deviations caused by the semi- 
classical approximation are less than 1% or This justifies the use of this approximation 
for cylindrical systems in a magnetic field. 

Figure 6 shows the results for a cylinder with iV = 200 sites along the circumference 
and M = 1000 sites along the cylindrical axis. The temperature is k^T = 10 -4 in units of t 
(around IK) and the density is 0.74 electrons per site. Since the polarisability is symmetric 
with respect to the point <fi = 1/2, as mentioned before, only the regime < <fi < 1/2 is 
shown. Obviously there is a huge effect of the magnetic field on the polarisability. 

These conclusions are based on a semi-classical approximation. In the next subsection 
we present full quantum calculations for the case where the electric field is perpendicular to 



9 



the cylindrical axis. Therefore, we postpone a complete discussion of the results to the next 
subsection where they do not suffer from any approximation in the context of our simplified 
model. 

B. Transversal electric field: Quantum calculation 

The Hamiltonian describing the system is now given by Eqs. fll2|) and ([14]). The advan- 
tage of this choice is to allow for a separation of variable: the motion of the electrons along 
the circumference and the cylindrical axis can be treated separately. Exact calculations can 
be done even for very long cylinders. Proceeding that way, the problem is reduced to a 
study of one-dimensional rings in applied electric and magnetic fields ||. 

We show in Fig. 7b the magnetic field dependence of the polarisability for a cylinder 
with iV = 101 and M = 100 and very few electrons on it, N e = 100 (this corresponds to 
a band filling of only 1%). In this example, the electric field is such that eaE = 10~ 3 t 
(linear regime) and the temperature is k{F = 10~ 5 t (T ~ 0.1K). We choose this example 
because it shows very clearly the main behaviours of the dielectric response. In particular 
we choose iV = 101 because this gives an illustration of the signature of the linear Stark 
effect. Together with this example the energies of the ground state and first excited state 
are shown without electric field as function of the magnetic flux (Fig. 7a). 

First of all the induced dipole moment as function of the magnetic field is periodic with a 
period of 1 and is symmetric with respect to 1/2; these symmetries are already apparent in 
the spectrum ((?]). Second, the induced dipole moment shows clearly two main characteristics: 
(i) we can notice the presence of small peaks at <fi = 0, 1/2, 1, (ii) the induced dipole moment 
is a discontinuous function showing several jumps separating plateau-like sections. Note that 
there is a slight curvature in the whole spectrum, which is related to the persistent current 
induced by the magnetic field. 

As noticed before, a magnetic field induces crossing between the energies of the ground 
state and the first excited state (figure 4 and 7a). As it can be seen in figures 7, there is a 
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one to one correspondence between level crossing at zero electric field and each kink in the 
induced dipole moment. 

At each crossing, the ground state of the system changes. These two states which are 
crossing respond differently to an applied electric field. Both produce a quadratic Stark 
effect but generally of different size. This explains why the induced dipole moment is not a 
continuous function of the magnetic field. 

Near the crossing points, the response of the system will depend on whether or not there 
is a nonzero matrix element of the electric field between the two states involved. If the 
matrix element vanishes, the picture described above is valid. On the contrary, if there is an 
interaction between those states, due to the degeneracy, the response will result in a linear 
(instead of quadratic) Stark effect with peaks of the induced dipole moment. More precisely, 
using the expression of the wave function @, the matrix elements of the dipole operator are 
calculated as 

< ^ m \d\^ p ^ q , >= eR5 pl)P±l 5 q ,, q (20) 

With this equation, it is easy to see that a linear Stark effect can occur only for a magnetic 
flux of = 0, 1/2, 1. These are precisely the values for which one obtains peaks in our first 
example (Fig. 7b). 

It is worth noticing that large enough Coulomb interactions can change drastically the 
selection rules (|20|). The same is true for distortions of a cylinder as it was shown in |8j for 
elliptical cylinders. Appearance of a linear Stark effect for different values of the magnetic 
flux could then give a way to quantify importance of electronic correlation effects or to 
characterise the geometry of the systems under consideration. 

At low enough temperature, one deals essentially with the spectroscopy of a few levels 
around the Fermi level. Therefore it is not surprising that the induced dipole moment for a 
system which is well ordered, well oriented and of well defined size (N, M), depends strongly 
on the electronic density. This is clearly seen from Figs. 8a and 8b. They are for the very 
same system considered before but for three different electron fillings N e = 100, 101 and 
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102. At such a low density, the relevant mean level spacing, AE, behaves like AE m 1/N 2 . 
Figure 8a is for a temperature lower by one order of magnitude than AE, while Fig. 8b 
is for a temperature higher by one order of magnitude than AE. With a typical value for 
t (t ~ 2eV) one can estimate a temperature of lO -1 ^ for the spectrum 8a and 10K for 
the spectrum 8b. With these figures we want to emphasise the unique sensitivity of the 
proposed measurements and also that it is the necessary to work at very small temperatures 
in order to extract maximum informations. 

The last examples in Fig. 9 is an illustration showing that the same kind of behaviour is 
expected in carbon nanotubes. It is for armchair (Fig. 9a) and zig-zag (Fig. 9b) nanotubes, 
when both are 1% away from half-filling and at T ~ 0.1K ||. 

Today it is possible to measure accurately very small variations in the real part of the 
dielectric function HQ- Therefore, the dramatic magnetic field effects on the static polar- 
isability of mesoscopic cylindrical systems, discussed here, could be measured and analysed. 
Several important informations about excited states could then be obtained. First, the po- 
sitions of the discontinuities and peaks should give informations about the energies of the 
excited states. The nature of the response - linear or quadratic Stark effect could be detected 
and therefore should give informations about the symmetry of the excited states. The mag- 
nitude of the response should also give informations about the coupling constants. Finally, 
the different curvatures observed in the whole spectrum could give informations related to 
the persistent current - or ring current - induced by the magnetic field. 

Until now we were discussing cylindrical systems, such as carbon nanotubes, with a 
diameter in the mesoscopic range. The diamagnetic current induced by the applied magnetic 
field is a ring current in a strict geometrical sense. In the following, we will discussed quite 
different systems closer to those discussed by Ehrenfest, Pauling and London [0-0 • Indeed, 
we consider finite graphene planes where the magnetic field induces a ring current running 
along the perimeter of these systems. It is very close to the earlier studies where ring currents 
in benzene or naphthalene molecules were considered. The mesoscopic graphene sheets are 
new materials of importance for which the magneto-polarisability [|l!| could give very useful 
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informations. 



IV. STATIC ELECTRIC POLARIS ABILITY UNDER MAGNETIC FIELD OF 

NANOGRAPHITE LAYERS 

Very recently a new type of magnetism similar to a spin glass state was found for acti- 
vated carbon fibres prepared at heating temperature over 1200 ° C ||20|| . These materials are 
supposed to be well described by disordered networks of nanographite layers with a charac- 
teristic length of some tens of Angstoms ||. It is worth noting that further increase of the 
heating temperature will drive the system towards complete graphitisation. 

On the other hand, recent theoretical works on graphite ribbons have shown the existence 
of edge states whose wave functions are localised in the near vicinity of the edges of the 
ribbons and whose energies form a band with a small dispersion centred in the middle of 
the gap pl| , |22| . Therefore, these states can have important physical consequences since they 
make a relatively large contribution to the density of state at the Fermi level. For instance, 
they are expected to affect novel electrical and magnetic properties such as a paramagnetic 



behaviour at low temperature |[22| . In the light of these theoretical works, it was suggested 



by Shibayama et al. that the edge states should be responsible for the spins observed in |20 
The spin glass behaviours could then be explained by invoking Heisenberg like interactions 
between these effective spins since the disordered network of nanographite layers gives rise 
to a random exchange interaction network. 

In this context, we want to stress in this section that measurements of the magneto- 
polarisability may be a very direct and efficient way to detect possible edge-states in activated 
carbon fibres or any other amorphous carbon material. This new way of probing electronic 
structure of carbon materials could help to support the hypothesis made in ||20|| . Moreover, 
we believe that this quantity could give a way to access characteristic length of graphitic 
regions in such compounds - quantity studied today using Raman spectroscopy p3|. 



For that purpose, we consider first rectangular finite planes of graphite (Fig. 10). These 
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simplest systems show two different kinds of edges: zig-zag and armchair edges. The one- 
band it electrons are again described by a nearest-neighbours tight binding Hamiltonian 



defined on the topologically equivalent brick- type lattice already used in [22 



H= J2 tat n,m a n',m' (21) 

<(n,m),(n' ,m')> 

where the summation is over nearest-neighbour sites located at (x, y) coordinates, (n, m) 
and (n', m'). 

As found in references pl| , p2| , our studies of ( plj) show the existence of edge states for 



finite rectangles of graphite. Moreover, in agreement with the results of pl| , |22[| for ribbons, 
the wave functions of the edge states are shown to be localised preferentially on the zig-zag 
edges. In particular the weight of these functions is negligible on the armchair edges. This 
remark could be of importance for real materials leading to the question about the stability 
of such states in layers with disordered shape. We postpone a discussion of this important 
question to a future publication. Finally, the energies of these states are near the middle of 
the gap; the closer it is to the middle, the more localised on the edges is the corresponding 
wave function. 

Next we apply a uniform magnetic field H perpendicular to the layer. In the tight binding 
model we then proceed once more by using the London substitution ||, t — > t( n ,m)->(n',m') = 

j 27r JL \' {n ' ' m>) dlA 

f e ch j( n ,m) ; where A is the vector potential. We write 

t(n,rn)— >(n+l,m) Lc ^ x ^ ' 

(22) 

ft \ I , -n = fp i T!/( n ' m ) 
b (n,m)^(n,m+l) vc - 

where , y x (n,m) and j y (n,m) are the appropriate line integral of the vector potential. Since 
the magnetic field is uniform, one must have (cf Fig. 10), again in units of 0o 

7„(n,m) - 7^(n + 2,m + 1)+ / x 

(23) 

j x (n, m + 1) + 7 x (n + 1, m + 1) - j x (n + 1, m) - j x (n, m) = 27r</>. 



For calculations we choose, as in the Landau gauge 
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7v(n,m) = 

(24) 

7a; (n, m) = n(f)m 

It should be noticed that the elementary magnetic flux is defined here with respect to the 
size of a unit cell, namely a benzene ring, which requires an extremely large field except for 
the edge states where the relevant unit is given by the size of the whole layer - corresponding 
to a much more reasonable magnetic field of some tens Tesla. In that sense, the edge-states 
define an effective rings. The idea is quite simple: only the edge states will react to small 
values of an applied magnetic field such that the corresponding characteristic length of the 
cyclotron radius is of the same order than the size of the layer. Subsequently, the signature 
of the edge states should be seen in the polarisability as function of magnetic field. 

We then add a static uniform electric field in the plane of the layer and calculate the 
magneto-polarisability for different electron-fillings. In this case, the dipole operator takes 
the following general expression d = e J2i ^i<4 a i where the site vector r$ is defined with 
respect to an arbitrary origin. We choose a system with a total of 76 benzene rings with 
10 along the zig-zag edges. We work at small temperature ksT = 10~Hij which could be 
estimated to T ~ 0.1K. 

The particular example studied here shows only two nearly degenerate edge states at 
zero energy. The edge states will be populated at half-filling but since the temperature is 
low compared to the energy gap, they are not occupied if one reduces the electron number by 
one. Figure (11) shows the magneto polarisability for these two cases. We see very clearly a 
quadratic dependence of the polarisability as function of the magnetic field for the half-filled 
case only. For the other case, no sensible effect are seen. We can then conclude that the 
edge states are responsible for this quadratic behaviour in accordance with the argument 
given above. 

To conclude this section, it should be stressed that the relative variations of the polaris- 
ability in Fig. 11 for the case with occupied edge states can be detected nowadays. Therefore, 
we believe it is possible to identify the existence of edge states in activated carbon fibres, 
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or other amorphous carbon materials, using measurements of the magneto-polarisability. 
Moreover, we believe that a careful analysis of the quadratic behaviour shown in Fig. 11 can 
help to gain important informations about these disordered materials such as the effective 
sizes of the graphite-like regions for instance. 



V. RING-CURRENT IN STRONGLY CORRELATED CYLINDRICAL 

NANOSTRUCTURES 

So far we didn't consider Coulomb interactions. Instead, we were assuming that they were 
included in the effective one-electron parameters of our tight-binding models. However, for 
quantitative results of the polarisation, it is needed to include explicitely those interactions. 
It is well known that the absolute value of the polarisation of organic compounds decreases 
strongly by including interaction corrections [p4|| . We will not attempt such studies here 
which will be the subject of a subsequent publication. Instead, we will consider several 
strongly interacting systems, rings or small cylinders, in the extreme limit of infinite on-site 
interactions, our purpose being to point out that Coulomb interactions can produce not 
only quantitative changes but, also, important qualitative changes in the behaviour of the 
polarisability. 



A. Fractional Aharonov-Bohm effect 



In the early nineties, several authors have studied the behaviour of persistent currents in 



mesoscopic Hubbard rings threaded by a magnetic field One of the most prominent 

result of these studies was obtained in the limit of infinite interaction, U — > oo. In this 
limit, Kusmartsev pointed out that, in the case of spin 1/2 fermions, the persistent current 
has a surprising short period given by A(f> = |p (in unit of <fio), where N e is the number of 



electrons ||25|| . He called this intriguing new phenomena, fractional Aharonov-Bohm effect. 
It is this result which we shortly describe here. 

The tight binding Hubbard model for a ring with N sites is of the form 
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H = tJ2(al+i,aan,ae i "' t> + h.c.) + h n ^h nA (25) 



n,a 



where n n>ff is the usual electron number operator. The magnetic flux is given by the 
expression (|2J). 

An advantage of the Hubbard ring is to be one of the strongly-correlated-electron system 
that is solvable by using the Bethe ansatz |27J. In the presence of a magnetic field the form 
of the wave function doesn't change 

-, Ne 

ty(x 1 ,x 2 ,...,x Ne ) = J2[Q,P)exp(iJ2 k Pj x Qj) ( 26 ) 
p i=i 

where P = (Pi, P 2 , P/vJ and Q = (Q\, Q 2 , Qn b ) are two permutations of (1, 2, N e ). 
The coefficients [Q, P] as well as (k\, k 2 , fcjvj are determined from Bethe's equations which 
in a magnetic field are modified by the addition of the flux phase <ft [ p5| - |26 . 



... , T n T-f .it sin kj — iX g — U/A, 

FV 3 YJ ^it sin kj - iXp + U/A' 



(27) 



^ iUmkj -iXp-U/A _ ^ iXg-iXp + U/2 

f}} it sin kj -iXp + U/A* ^iX a - iX -U/2 } [ ' 

N e is the number of charges i.e. the number of electrons, and N s is the number of down- 
spins. The eigenstates of a ring are characterised by the momenta kj of the charges and the 
rapidities A/3 of spin waves. The energy of the system is given by 

E = 2tJ2 cos kj (29) 
i=i 

In the limit of infinite interaction U — > 00 Bethe's equations are greatly simplified to 
give the following results 

2tt 1 Ns 

^ = ¥ (^+0+TF-£ J a) (30) 

e a=l 

where the quantum numbers Ij and J a are connected with charge and spin degrees of freedom 
respectively. They are either integers or half odd integers, depending on the parities of the 
numbers of down- and up-spin electrons: 
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lj = -—{moal), J a = {modi). (61) 

The equation (|30D is identical to the one describing a set of non-interacting spinless 
fermions on a ring threaded by a flux <f> + ^- J2a=i J a- The second term of this expression 
is the effective magnetic flux created by the spin excitations. It comes from the fact that 
there is an additional phase shift between the momenta of the different electrons due to the 
repulsive interaction. Because of this shift, the periodicity of the Aharonov-Bohm flux takes 
a fractional value at the infinite interaction limit. It is easy to check this statement by using 
Eqs. (EUp and (BDp and by appropriately choosing the quantum numbers (PTI) in order to 



minimise the energy of the system. One example is shown in Fig. 12a for a ring with 8 sites 
and 2 electrons, one up and one down. Another example is shown in Fig. 12b for the very 
same ring but for 5 electrons, 3 up and 2 down. The fractional Aharonov-Bohm effect is 
clearly seen in these examples. 

It is important to stress once more that the spin plays an essential role in this phenomena. 
This can be seen in (|30|) where the quantum numbers associated to the spin degree of freedom 
appear as responsible of the fractional AB effect. Without this term we recover the known 
result for free fermion rings (^). In particular, for spinless fermions, the periodicity with the 
flux is 1. 

In the case of infinite interaction or spinless fermions, we can use the following argument 
to obtain the correct flux periodicity and which stresses the important role played by the 
spin. Let us take a particular configuration with M electrons, \x%, x%, xm >■ One may 
ask the question of how many hops of the particles in the same direction should take place 
in order to recover the starting configuration. In mathematical terms we may write 

(e^TO^e^^^^ (32) 

where a is a phase term cause by the Pauli's principle. Xj, are general coordinates of the 
electron i taking account of the spin when necessary. Tj are the operator of translation 
by one site in the clockwise direction (for instance) of the electron i. iVj is the number of 



one-site-hopping of the particle i. The answer to this question depends on whether or not 
we consider the spin degree of freedom. 

In the cases we are interested in, the fermions are locked together and form a kind of 
dynamical Wigner crystal. By moving on the lattice, an electron will push the others in the 
same direction. For the case of spinless fermions, since all particles are indistinguishable, it 
is easy to see that N± + N2 + ... + Nm — N to get back to the initial configuration. The 
same is true for an assembly of free spin 1/2 fermions. On the other hand, one come back 
to the same configuration by interchanging the particles which, due to the Pauli's principle, 
give the following statistical flux: a = —e mM . The immediate conclusions are, first, that 
the periodicity of the energy with the flux will be 1 (as for the free electron case) and, 
second, there will be a parity effect which discriminates between the cases with even and 
odd number of particles. Indeed, it was shown by Kusmartsev [[25|] (and references therein), 
that a system of interacting spinless fermions on a ring behaves as a paramagnet for even 
number of particles and as a diamagnet for an odd number of particles. This property is 
sometimes known as Leggett's conjecture 



For the case of spin 1/2 fermions, the conclusions are completely different. Indeed, one 
easily sees that Ni + N 2 + ... + Nm = N.M and a = 1. Each particle should go round the 
ring and therefore, the statistical flux is just unity. Consequently, the effective charge of 
the system will be Me and the periodicity of the energy will be l/M - this is the fractional 
Aharonov-Bohm effect discovered by Kusmartsev [p5|j . Moreover, there is no parity effect: 
indeed, it is known that the system is always diamagnetic for infinite interactions |29| . 



On the basis of these analysis, let us now consider strongly-interacting cylindrical systems 
in some limiting cases where conclusions could be reach following the same argumentation. 

B. Strongly interacting rings with one electron per ring 

Let us first consider a stack of isolated rings with only one electron per ring. The rings 
are interacting following the model 
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H = t^2(a ] n+ltm a n ^ m e 1 ^ ^ + h.c.) + U ^ i,mi,m+i (33) 



n,m 



where there is a nearest-neighbour interaction between two electrons when they are on top 
of each other. This Hamiltonian could help to describe stack of large aromatic molecules 
with only one electron in each in the LUMO. It could also be an effective model to describe 
interacting edge-states in stack of finite graphene planes. This model was already used in 
[|T0f| in the continuum limit, where an interesting new parity effect was pointed out to occur 
in the dielectric response of one-dimensional stack of quantum rings. Here, we are only 
interested by the periodicity of the ground state energy with the magnetic flux in the limit 
of infinite interaction. 

Once again, starting from one particular electronic configuration, we are asking the 
same question about the number of particle hopping in the same direction one should do 
to come back to it. It is then possible to write down a similar expression than ( [32]) with 
Ni + N 2 + ... + Nm = N.M and a = 1 where M is the number of rings and N, the number 
of sites in each ring. Consequently, one can conclude that the periodicity is 1/M and that 
there is no parity effect. The label of the ring acts as a spin degree of freedom. It is worth 



mentioning that it exists a rigorous mathematical proof of these properties 3C 



C. Strongly interacting cylinders 

Last, we are considering a model of interacting cylinders 

H t'Yln,m,a(fln+l,rn,(jQ'n,m,<j& N ^ ~\~ h.C.) + t X/n,m,o-( a n,m+l,cr a n,m,(T 4" /l.C.) + 

(34) 

with again U — > oo. Again we are asking the same question about the periodicity formulating 
it in terms of equation (p2|). One has to distinguish between the case of spin 1/2 and spinless 
fermions. 

• For spinless fermions, for a non vanishing t', one has N% + N 2 + ... + Nm = N and 
a = 1; the periodicity is one and there is no parity effect. 
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• For spin 1/2 fermions, whatever is the value of t', one has N\ + N2 + ... + Nm — N.M 
and a — 1; the periodicity is 1/M and there is no parity effect. 

In conclusion, there is an abrupt transition from a fractional Aharonov-Bohm state to 
an usual Aharonov-Bohm state going from t' = to t' 7^ in the case of spinless fermions; 
on the contrary, there is no such transition for spin 1/2 fermions. 

VI. CONCLUSION 

The main purpose of this review is to point out the strong magnetic field dependence of 
the static polarisability of nanoscopic cylindrical organic materials when the magnetic field 
is parallel to the tube axis. This was done by studying two kinds of system. 

On one hand we have considered real cylinders in the form of rolled square lattices but 
also, rolled honeycomb lattices, more specifically, achiral carbon nanotubes. They were 
studied with both, semiclassical approximation [7|] and exact diagonalisation ||. For all 
these cases, the polarisability was shown to present very complex structures as function of 
the magnetic field in which one can identify two different characteristics (cf Fig. 7): (i) 
the polarisability is a non-continuous function with sudden jumps separating plateau-like 
regions (ii) additionally, small peaks may appear for special values of the magnetic field in 
place of jumps. A full understanding of these complicated behaviours was given by following 
the behaviour of the ground state by increasing the magnetic field: due to the Aharonov- 
Bohm effect, many changes of ground state occur and for some values of the magnetic field 
accidental degeneracies happen where the ground state becomes two fold degenerate. A one 
to one correspondence is found between the accidents in the polarisability and the accidental 
degeneracies of the ground state. Each plateau-like region of the polarisability corresponds 
to a quadratic Stark effect with a coefficient proper to the corresponding magnetic-field 
induced ground state. Each peak corresponds to a linear Stark effect appearing at accidental 
degeneracies when there is direct coupling between the two states involved (cf Fig. 7). 

On the other hand we have studied finite plane of graphene which support edge states. 
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The wave functions of these states have the peculiarity to be mostly localised on the perime- 
ter of the plane defining kind of effective rings. The presence of these states in real materials 
is a important question since they are able to influence new physical properties [2(| and, 
we have shown that it is possible to use measurements of the static polarisability in an 
applied perpendicular magnetic field to detect possible edge states in amorphous carbon 
materials such as carbon fibres. Moreover we believe that a careful analysis of the quadratic 
behaviour we have pointed out in section IV (Fig. 11), should give a way to determine the 
characteristic length of the graphitic region in these disordered systems. 

In a last part of this review, we have briefly described results of some strongly interacting 
models in the limit of infinite interactions. The main purpose of this section was to pointed 
out that the interaction can result in new phenomena such as the fractional Aharonov-Bohm 
effect which appears for Hubbard rings with spin 1/2 fermions but also in stacks of uncouple 
rings with one electron per ring. An explicit treatment of the electron-electron interactions 
is therefore needed and, this is now our main direction of investigation. 
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FIG. 2. Cylinder made of rolled square lattice placed in a uniform magnetic fields, H, parallel 
to the cylindrical axis. A uniform electric-field is applied parallel to the cylindrical axis (section 
Ill-a) or perpendicular to it (section Ill-b) 
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FIG. 3. Lowest eigenvalues (Eq. (0)) for a cylinder with N = 100 sites along the circumference 
and M = 100 sites along the axis, as function of flux eft through the cylinder (in arbitrary units). 
Level crossing appear at values of <j) depending on the cylinder geometry. 
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FIG. 4. Energies of the ground state (full line) and of the first excited state (dotted line) as 
function of the magnetic flux, <p, for a cylinder with N = 100, M = 100 and N e = 1000 (the 
references are the energies without magnetic field). 
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FIG. 5. Induced dipole moment as a function of electrons per site for a chain of N = 201 



sites calculated with the exact quantum mechanical expression ([15]) (dashed lines) and with the 
semi-classical expression ( |19| ) using the Thomas- Fermi approximation (solid lines), (a) and (b) 
correspond to temperatures ksT = 0.01 and 0.05, respectively (in units of the transfer integral). 
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FIG. 6. Dielectric response (e'(H) — e'(0))/e'(0) for a model square-lattice system with N = 200, 
M = 1000 in an axial magnetic field. The temperature is ksT = 10~ 4 (in units of the transfer 
integral), and a = 1AA. The electron density is 0.74 per site. 
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FIG. 7. Cylinder with N = 101, M = 100 and iV e = 100, the numbers of site along the 
circumference, along the cylindrical axis and number of electrons (respectively). (a) Energies of the 
ground state (full line) and of the first excited state (dotted line) as function of the magnetic flux, 
<p, without electric field (the references are the energies without magnetic field), (b) Polarisability 
(P(T, 0) - P(T, 4>))/P(T, 0), at k B T = 10" 5 t and for v = eER = W~ 3 t as function of the magnetic 
flux; t is the hopping integral defined in (0), E the electric field and R the radius of the cylinder. 
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FIG. 8. Polarisability of a cylinder with N = 101 and M = 100 as function of the magnetic flux 
(j) in arbitrary units (a.u.). (a) For N e = 100, 101, 102 from the bottom to the top at a temperature 
of k B T ~ lO" 1 ,^. (b) The same but for k B T ~ 10K. 



32 




33 




(n,m+l) 




o 
cr 

&■ 
>-t' 

ffl 

D. 

ro 







(n+l,m) 










(n,m) 



















FIG. 10. (a) Example of a nanographite plane with rectangular shape. The boundaries have two 
different kinds of edges: the armchair and the zig-zag edges. The wave function of the edge states 
are preferentially localised on the zig-zag edges, (b) Corresponding Brick-type lattice equivalent 
to the honeycomb lattice shown in (a). The sites are labelled with the of (n,m) indices 
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FIG. 11. Relative variation of the polarisability, (P(T, <p) — P(T, 0))/P(T, 0), of a nanographite 
plane with 76 benzene rings as a function of the magnetic flux <f>, in units of (j>o, for T ~ 0.1K. 
4> = 10 -5 corresponds to a magnetic field of IT approximately. The full line is obtained for the case 
of half-filling; the dotted lines is for the case where we have reduced the electron number by one 
with respect to half-filling. The edges-states are populated only at half-filling and are responsible 
for the quadratic dependence seen in this figure. 
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FIG. 12. Ground state energy as function of the magnetic flux <p for Hubbard rings with infinite 
interaction for a ring with 8 sites and (a) 1 electron up and 1 electron down (b) 3 electrons up and 
2 down. In both cases the fractional- Aharonov-Bohm effect is clearly seen. 
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